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Abstract 


An accurate accounting of the surface emissivity is 
important both in the retrieval of surface temperatures and in 
the calculation of the longwave surface energy budgets which 
are derived from data collected by remote sensing instruments 
aboard aircraft and satellites. To date, however , high quality 
surface emissivity data have not been readily available for 
global applications. As a result, many remote sensing and 
climate modeling efforts have assumed the surface to radiate 
as a blackbody (surface emissivity of unity). 

Recent measurements of spectral reflectances of surface 
materials have clearly demonstrated that surface emissivities 
deviate considerably from unity , both spectrally and 
integrated over the broadband. Thus , assuming that a surface 
radiates like a blackbody can lead to potentially significant 
errors in surface temperature retrievals in longwave surface 
energy budgets and in climate studies. Taking into 
consideration some recent spectral reflectance measurements , 
we have constructed global maps of spectral and broadband 
emissivities that are dependent on the scene (or surface ) type. 
To accomplish our goal of creating a surface emissivity map , 
we divided the Earth* s surface into a 10' lot. X 10' Ion. grid, 
and categorized the land surface into 18 scene types. The first 
17 scene types correspond directly to those defined in the 
International Geosphere Biosphere Programme (IGBP) 
surface classification system. Scene type 18 has been added to 
represent a tundra-like surface which was not included in the 
IGBP system. Laboratory measurements of the spectral 
reflectances for different mineral and vegetation types were 
then associated, individually or in combination, with each of 
the 18 surface types, and used to estimate the emissivities of 
those surface types. Surface emissivity maps were generated 
from the band-averaged laboratory data for 12 longwave 
spectral bands (> 4.5 pm) used in a radiative transfer code as 
well as for the NASA's Clouds and the Earth's Radiant Energy 
System (CERES) window channel band (8-12 pm). The 
spectral emissivities for the 12 spectral bands were 
subsequently weighted using the Planck function energy 
distribution to calculate a broadband longwave (5-100 pm) 
emissivity. The resulting broadband emissivities were used 
with a surface longwave model to examine the differences 
resulting from the use of the emissivity maps and the 
blackbody assumption. 



Introduction 


Measuring the longwave (LW) radiation budget at the Earth’s surface is a critical part of 
NASA’s Clouds and the Earth’s Radiant Energy System (CERES) project (Wielicki et al. 

1996) . Such LW measurements will foster a better understanding of the energetics of the 
Earth’s atmosphere- surface system. Accurately characterizing surface emissivity (e) is 
important for correctly determining the longwave radiation leaving the surface and for 
retrieving surface temperature from remote sensing measurements (Wan and Dozier 1996; 
Kahle and Alley 1992; Kealy and Hook 1993). Frequently, past studies have assumed the 
emissivity to be unity when determining surface temperature and longwave emission. To 
improve the accuracy of retrieved satellite products, global maps have been created of surface 
emissivities for 12 longwave spectral bands of the Fu-Liou radiative transfer model (Fu and 
Liou 1992), for the CERES window channel (8-12 jam), and for the broadband longwave 
region. 

Ideally, field measurements of a wide range of surface types (e.g., soils, crops, forests, 
grasslands, and semi-arid) would be sufficient to construct surface emissivity maps. Accurate 
in-situ measurements of emissivity, however, are very difficult to obtain because the 
parameters which influence apparent emissivity, namely the surface temperature and 
atmospheric state, are highly variable quantities and are difficult to measure. In addition, the 
spatial coverage of the available field measurements is insufficient for global studies. Remote 
sensing measurements from satellites could be used to retrieve surface emissivities, but that 
requires concurrent temperature measurements on the ground as well as detailed knowledge 
of atmospheric absorption and scattering. With a single measurement of surface temperature 
or emissivity, the problem is undetermined. Many methods have been used to approach this 
problem. There has been some success (Van de Griend and Owe 1993; Olioso 1995) in 
relating the thermal emissivity to the Normalized Difference Vegetation Index (NDVI) from 
the Advanced Very High Resolution Radiometer (AVHRR). Several algorithms have been 
developed for use with the new generation of satellite instruments. The Advanced 
Spacebome Thermal Emission and Reflection Radiometer (ASTER), the Moderate Resolution 
Imaging Spectroradiometer (MODIS), and the Process Research by Imaging Space Mission 
(PRISM) instruments will all use temperature-emissivity separation algorithms to retrieve 
surface temperature and surface emissivities from remote measurements (Caselles et al. 

1997) . The new generation of satellite instruments and the current thrust of field experiments 
will help to fill the gap in our knowledge of surface emissivity. 

Algorithms developed for the CERES processing use surface emissivity to determine the 
longwave radiation budget at the Earth’s surface. As a consequence, surface emissivity maps 
were needed as soon as the CERES instrument began taking measurements. Nevertheless, 
surface emissivities from other EOS instruments will not be made available until 3 to 5 years 
after the launch of the first CERES instrument. We have, therefore, created surface emissivity 
maps from laboratory measurements. Such maps will constitute a viable source for surface 
emissivity data until the EOS surface emissivity measurements become available to create 
more advanced maps. 

The emissivity of the ocean surface is known to vary with the viewing zenith angle and the sea 
state. When determining the sea surface temperature from space, this variation is important. 
The effects of viewing angle and sea state on emissivity have been both modeled and 
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measured (Masuda et al 1988; Smith et al 1996; Wu and Smith 1997); however, it has been 
found that at near-nadir viewing angles both spectral and broadband surface emissivities are 
nearly constant with respect to sea state. Since the present global emissivity maps have been 
developed for nadir viewing conditions, we do not need to consider the effect of sea state on 
emissivity. More advanced versions of the surface emissivity map will consider sea state 
whenever zenith angle effects are included. 

To accomplish the goals of the CERES project, the Surface and Atmospheric Radiation 
Budget (SARB) group in the CERES experiment created a global scene type map on a 10 
grid (Rutan and Charlock 1997). The scene types for this map were adopted from the 
International Geosphere Biosphere Programme (IGBP) classification for the Earth’s surface 
with an additional scene type for tundra. Global maps of broadband and spectral albedos 
necessary for calculating the shortwave (SW) portion of the radiation budget were developed 
based on the IGBP scene types. The present emissivity maps have, for the sake of 
consistency, utilized the same IGBP scene types and thus are compatible with the albedo 
maps. The MOD1S instrument team has also proposed using a classification-based method to 
determine surface emissivity (Snyder et al. 1998). For instance, the MOD1S team will use 
surface emissivity data for MOD1S channels 31 and 32 in order to retrieve land surface 
temperature. Laboratory measurements of the spectral reflectances of several different 
vegetation and soil types found in the Johns Hopkins Spectral Library (Salisbury and D’Aria 
1992a) have been used to compute spectral emissivities for each of the 18 surface types. The 
spectral reflectance data have been used to calculate band-average emissivity in each of the 12 
spectral bands used in the Fu-Liou radiative transfer model. The 12 band-average emissivities 
were combined into broadband emissivity by weighting with the Planck function energy 
distribution. The CERES window channel has also been subdivided into three bands. The 
emissivities of the three bands were then combined into a window emissivity in the same 
manner as the broadband. 

The global emissivity maps presented herein are the first to handle surface emissivity on a 
scene dependent basis. Until surface emissivity measurements become available for a 
representative fraction of the earth, it is beneficial to have global maps of surface emissivity 
based on surface type using presently available measurements. The global surface emissivity 
maps will be modified as more laboratory and field data become available. In addition, data 
received from the MODIS and ASTER instruments will also be incorporated. Even after we 
make available more advanced versions of the emissivity maps which consider zenith angle 
effects, seasonal effects, and wider variety of surface types, the present emissivity map should 
still prove useful for use by models where more complex maps are not warranted. 

Theory and Background 

Laboratory measurements have been made for the spectral reflectances of a wide variety of 
surface materials (Salisbury and D’Aria 1992a). Since the CERES processing algorithms 
require surface emissivities, the measured reflectances were transformed into emissivities by 
applying energy conservation and Kirchhoff s Law. 

For a process involving absorption, reflection, and transmission, the total energy, normalized 
to unity, is partitioned as: 
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Ax + Rx + T x - 1 . 


(1) 


The present study assumes the transmittance {T ; ) of the surface to be zero. Applying 
Kirchhoff s Law which states that the absorptance (A a ) is equal to the emittance (e A ) under 
conditions of thermodynamic equilibrium yields a straightforward relationship between 
reflectance (R A .) and emissivity: 

= 1 - Ra- (2) 

Figure 1 shows emissivity as a function of wavelength derived from laboratory measurements 
for nine materials representative of a large percentage of the Earth’s surface. 





Figure 1. Emissivity from Laboratory Measurements for 9 Materials. 
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Fu and Liou (1992) developed the radiative transfer code which is being used by the SARB 
group in the CERES experiment. The code uses a del ta-four- stream (Liou et al. 1988) 
approach for scattering and a set of correlated k-distributions for absorption in 18 distinct 
spectral bands. This code divides the shortwave region (0.2-4. 5 jam) into six bands, and the 
longwave (>4.5 jam) into 12 bands. The longwave bands, which cover the spectral range of 
interest in this report, are presented in Table 1. 


Table 1. Specification of Fu-Liou bands in the longwave region. 


Band 

Wavelength (jam) 

Wavenumber (cm 1 ) 

1 

4.5-5. 3 

2200-1900 

2 

5. 3-5. 9 

1900-1700 

3 

5.9-7. 1 

1700-1400 

4 

7. 1-8.0 

1400-1250 

5 

8.0-9. 1 

1250-1100 

6 

9.1-10.2 

1100-980 

7 

10.2-12.5 

980-800 

8 

12.5-14.9 

800-670 

9 

14.9-18.5 

670-540 

10 

18.5-25.0 

540^100 

11 

25.0-35.7 

400-280 

12 

>35.7 

280-0 


The choice of wavelength ranges for the bands in a radiative transfer model is based 
primarily on the absorption/emission characteristics of the atmospheric constituents. 
However, radiative characteristics of the surface in these bands are just as important for 
analyzing the radiation measurements. A brief description of the absorption/emission 
processes in each of the Fu- Liou bands, and the resulting relationship betwen emissivity and 
outgoing LW radiation are presented below. 

Fu-Liou Bands 

The first band (4. 5-5. 3 jam) is fairly transparent; however, the energy emitted from the 
surface in this spectral range is quite small. As a result, this spectral range does not provide a 
great deal of information concerning surface and lower tropospheric properties. Bands 2 and 
3 (5. 3-5. 9 jam and 5.9-7. 1 jam) encompass the very strong 6.3 jam band of water vapor. 
Thus, very little energy emitted from the surface is transmitted through the atmosphere in this 
wavelength range. Band 4 (7. 1-8.0 jam) involves the important absorption due to the minor 
trace gases CH 4 and N 2 0, as well as additional water vapor absorption. The CH 4 absorption 
band centered at 7.7 jam, and the N 2 0 absorption band centered at 7.8 jam are important 
contributors to the atmospheric greenhouse effect. Since the 7. 1-8.0 jam spectral range 
allows some emitted surface radiation to escape to space, the surface emissivity becomes 
important to TOA measurements in this spectral range. When taken together, bands 5, 6 and 
7 (8.0-9. 1 jam, 9.1-10.2 jam, and 10.2-12.5 jam) correspond very closely to the CERES 
window channel. Since the rationale behind subdividing this spectral range is valid both for 
running the Fu-Liou model and for processing the CERES measurements, these intervals will 
be discussed in detail when we discuss the CERES window channel. Note, the spectral regions 
are by far the most transparent in the infrared. Thus, outgoing radiances in these regions 
coressponding to the CERES window channel are strongly affected by the surface emissivity. 
Bands 8 and 9 (12.5-14.9 jam and 14.9-18.5 jam) are characterized by the very strong 
15 jam band of C0 2 . While there is also a small amount of water vapor absorption within 
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these spectral intervals, it is the C0 2 absorption which dictates the widths of these spectral 
bands. Absorption and emission due to C0 2 in the 15 pm band provides an important 
source of radiative cooling throughout the atmosphere. The spectral interval from 12.5-14.9 
pm allows some radiation to escape to space, most notably near the short wavelength end. 
The spectral interval from 14.9-18.5 pm is essentially opaque to surface emission. Bands 10, 
11 and 12 (>18.5 pm) represent the water vapor pure rotation band. Because the pure 
rotation band of water vapor covers a rather wide spectral range and the absorption is far 
from uniform in nature, this spectral range has been subdivided into the three intervals. For 
conditions where water vapor burden is very low (< 3 kg m" ), the 18.5-25.0 pm band 
becomes relatively transparent, and is sometimes referred to as the “dirty window” or the 
“polar window” band. Thus, surface emissivity may be important in the 18.5-25.0 pm 
band. 

CERES Window channel 

The CERES window channel measures the thermal infrared energy emitted from the Earth 
within the spectral range between 8 and 12 pm. The CERES window channel is nearly ideal 
for measuring the energy emanating from the proximity of the Earth’s surface because its 
wavelength range corresponds to the most transparent part of the infrared spectrum. These 
measurements provide valuable information concerning atmosphere-surface interactions. 
Recent measurements (Salisbury and D’Aria 1992a) have demonstrated that the emissivities 
of typical terrestrial materials (soils, vegetation, water, etc.) can be significantly less than unity 
and are quite variable over the spectral range of the CERES window channel. Moreover, 
while the window region is relatively transparent, within this spectral range there exists a 
significant amount of highly nonuniform molecular absorption due to F^O and 0 3 , as well 
as a host of minor trace species (Kratz and Rose 1999). Furthermore, the opacity of thin 
cirrus clouds tends to increase significantly toward both longer and shorter wavelengths 
within this spectral range (see e.g., Figure 2 of Prabhakara et al 1993). The distribution of 
these processes within the CERES window channel has prompted the subdivision of the 
CERES window channel into three distinct subintervals, as shown in Table 2. Data from these 
3 spectral intervals will ensure accurate modeling of top-of-atmosphere (TOA) radiation in 
the CERES window channel. 


Table 2. Specification of the sub-intervals of the CERES window channel. 


Band 

Wavelength (pm) 

Wavenumber (cm 1 ) 

1 

8.0-9. 1 

1250-1100 

2 

9.1-10.2 

1100-980 

3 

10.2-12.0 

980-835 


The subinterval covering the 8.0-9. 1 pm range is associated with the very strong asymmetric 
stretching fundamental bands (reststrahlen bands) of quartz. Within this subinterval, bare 
soils have characteristically low emissivities. In contrast, the highest emissivities for senescent 
leaves occur within the 8.0-9. 1 pm spectral range. Moderately weak absorption due to F^O, 
as well as minor absorption due to O3 N 2 0, and CH 4 also occur within the 8.0-9. 1 pm 
spectral range. In addition, thin cirrus clouds possess a relatively high opacity in this 
subinterval. For the subinterval covering the 9.1-10.2 pm range, upwelling TOA flux is 
strongly affected by the 9.6 pm band of 0 3 . The presence of this strong 0 3 band along with 
some additional absorption due to H 2 0 and C0 2 prevents direct sensing of the surface by 
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satellites within the 9.1-10.2 |Lim range. Thin cirrus clouds possess a somewhat lower opacity 
in this subinterval as compared with the other two CERES window channel subintervals. The 
last subinterval, (10.2-12.0 pm), has surface emissivities for bare soils that tend to be 
relatively high. In addition, the highest emissivities for green foliage occur within this 
spectral range. The moderately weak H 2 0 absorption, as well as the smaller, yet significant 
contribution from the 10.4 pm band of C0 2 also affect this interval. The opacity of the 
cirrus clouds in this subinterval is comparable to that of the first subinterval. 


Method of Band Averaging and Weighting 

The emitted spectral radiance L x at wavelength X from a surface at temperature T s is 
calculated by multiplying the Planck function, B X (T S ) by the spectral emissivity e x 

h(T s ) = e x B x (T s ). (3) 

The Planck function, B x (T s ) represents the radiance emitted by a blackbody at a wavelength 
X and a surface temperature T s . 

Integrated over all wavelengths: 


n\B l {T s )dX = oT* , (4) 

0 

where a is the Stephen-Boltzman Constant (5.67 X 10 8 W m’ 2 K’ 4 ). 

The band-average emissivity in band i, from Tables 1 or 2, is defined by: 


ha. 


<i, upper) 


je x B x (T)dX 


— A a i 

£; = 


^ (i, upper) 

\B x {T)dX 


(5) 


The Planck function term in equation (5) can be taken out of the weighting process without 
introducing significant error. Such a strategy can be used for the following reasons. First, 
the wavelength dependence of the Planck function is relatively weak for the small intervals 
considered by Table 1. Second, the temperature dependence of the emissivity is usually very 
small for most surface materials. Even for coarse sands, which show the most variation, the 
band-averaged emissivity in the 3.5-4.25 pm band changes only 0.004 as the temperature 
changes from 240 to 320 K (Wan and Dozier 1996). Third, the spectral dependence of the 
surface emissivity is not generally correlated with the spectral dependence of the Planck 
function. Thus, the band emissivity calculated from laboratory reflectance spectra of a pure 
material becomes 
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(6) 


X 

^ (t, upper ) 

je x ca 

X 

— (i, lower ) 

~ 3 

^ (i, upper ) 

J dX 

^ (i, lower) 

which is taken to be independent of surface temperature because, as noted previously, the 
temperature dependence is usually very small. Equation (6) is used to calculate 

band-averaged emissivity. Figure 2 shows the location of the Fu-Liou spectral bands relative 
to the laboratory measured emissivity of a conifer sample. 



Figure 2. Locations of Fu-Liou Bands from Table 1 Relative to Laboratory Measurements of Conifer 
Sample. 

In general, the laboratory measurements spanned the wavelength range of 2-16 jam. At 
wavelengths greater than 16 pm, the emissivity was extrapolated, i.e., the measured emissivity 
in the interval closest to 16 pm was replicated to fill the remaining bands where data were not 
available. 

Using equation (6), we calculated 12 band-averaged emissivities for the Fu-Liou bands and 3 
band-averaged values for the CERES window channel. To calculate a broadband emissivity 
from the band-average emissivities, the Planck function was used to energy weight each of the 
12 band-average emissivities. The weighted values were then combined into a broadband 
emissivity. A weighting factor (wj) was calculated for each of the 12 bands as 
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(7) 


® (i, upper) 


B m d(0 


wfi = 


® (i, lower) 

“2200 


B C0 dco 


where, B C0 dco -B x dX , w/- is the weighting factor for a band and © Upper and © Uower are the upper 
and lower wavenumbers for each spectral band from Table 1. It is more convenient to use 
wavenumber space in the infrared to integrate over small intervals of the order of lOcmfl. 
Therefore, the Planck function was expressed in terms of wavenumber and integrated. The 
conversion factor from wavelength (pm) to wavenumber (cm’ ) is: ©(cm’ )=10000/A, (pm) 
[e.g., 1 pm= 10000 cm’ ]. In this manner a weighting factor was calculated for each band. 
The broadband emissivity, e^was then calculated by 


i = 12 

(8) 

i= 1 

Differences in temperature had no significant effect in the broadband calculation except for 
some minor effect in the case of quartz sand; therefore, a temperature of 288 K, which is 
representative of the average surface temperature of the Earth, was chosen for the calculation 
of the weighting factors. Note that a variation in temperature from 263 to 313 K resulted in a 
change of 0.011 in broadband emissivity for quartz. The change in emissivity of vegetation 
was 0.002 for the same variation in temperature. The weighting factors described in equation 
(7) were also used to combine the emissivities of the sub-intervals of the CERES window 
channel emissivity. 

Data 

Until recently, no emissivity data for vegetation were available in the thermal infrared region. 
Because of improvements in detector technology and measurement methods, measurements 
of spectral reflectivity have become available for various land cover types. We used the data 
from the Johns Hopkins Spectral Library. The ASTER team is also using this spectral library 
and has created an easily accessible database located at: 

http: // specli b. j pi .nasa. gov/ 

Data consist of laboratory measurements of reflectance of various types of vegetation, soils 
and snow and ice. 

Information on the measurement techniques is available at: 

h ttp : //spec li b . j pi . nasa . gov/ document s/j hu_de sc . htm 

http://speclib.jpl.nasa.gov/archive/JHU/becknic/vegetation/vegetation.txt 

http://speclib.jpl.nasa.gov/archive/JHU/becknic/water/snow&ice.txt 
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All spectra in the Johns Hopkins Spectral Library were measured under the direction of John 
W. Salisbury. Two similar instruments were used to record reflectance in the infrared range 
between 2.08-15 jam. Both are Nicolet FTIR spectrophotometers and both have a 
reproducibility and absolute accuracy better than plus or minus 1 percent over most of the 
spectral range. The data were quality checked at Johns Hopkins University (JHU). The 
instruments record spectral data in wavenumber space where both wavenumber accuracy and 
spectral resolution are given in wavenumbers (cm’ ). Wavenumber accuracy was limited by 
the spectral resolution, which yields a data point every 2 wavenumbers (cm’ ) for these 
measurements. The x-axis was changed from wavenumbers (cm’ ) to wavelength (jam) for all 
of the data before use in the present calculations. 

Spectra of vegetative canopy are not readily measurable. The lack of availability of field 
spectrometers and the effect of atmospheric absorption create difficulties in making field 
measurements. These problems are being overcome (Snyder and Wan 1996) and the results 
from the relevant studies will be included in future releases of the emissivity maps. 
Laboratory measurements have been made of simulated canopies from which the present data 
were derived. 

Measurement of the spectra of many different types of vegetation showed that conifer 
needles, deciduous tree leaves, and grass blades all have a very low reflectance (high 
emissivity) throughout the thermal infrared range. Because of the low reflectance and small 
spectral variation, one typical deciduous leaf spectrum was chosen to represent all deciduous 
species, one conifer to represent all conifers, and one grass species to represent all grasses. In 
a canopy, Salisbury and D’Aria (1992a) have noted several factors which combine to result in 
an emissivity quite close to unity. 

Assignment of surface types 


The surface types used are those from the IGBP from Bel ward and Loveland (1996) with the 
addition of tundra as type 18. Table 3 presents the 18 surface types. 

Table 3. International Geosphere Biosphere Programme Global Land Cover types. 


Type ID 

IGBP Type 

1 

Evergreen Needleleaf Forest 

2 

Evergreen Broadleaf Forest 

3 

Deciduous Needleleaf Forest 

4 

Deciduous Broadleaf Forest 

5 

Mixed Forest 

6 

Closed Shrublands 

7 

Open Shrubland 

8 

Woody Savannas 

9 

Savannas 

10 

Grasslands 

11 

Permanent Wetlands 

12 

Croplands 

13 

Urban 

14 

Cropland/Mosaic 

15 

Snow and Ice 

16 

Barren 

17 

Water Bodies 

18 

Tundra 
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A more complete description of the 18 surface types is presented in Appendix A, which has 
been adopted from Table 1 in Belward and Loveland (1996). 

Figure 3 illustrates the global distribution of the 18 surface types given by Table 3. The data 
is presented on the 10' grid used by CERES. 

As noted previously, the JHU spectral library contains reflectance measurements of a variety 
of surface materials. For the purposes of this study, 10 surface materials, individually or in 
combination are taken to be representative of the 18 surface types. The 10 surface materials 
used from the spectral library are: grass, conifer, deciduous, fine snow, medium snow, coarse 
snow, frost, ice, seawater, and quartz sand. The emissivity spectra of nine types of surface 
cover are shown in Figure 1. Table 4 shows how the 10 surface materials are associated with 
the 18 surface types. 


Table 4. Assignment of laboratory measurements to surface types. 


Type ID 

IGBP type 

Spectral library 

1 

Evergreen Needleleaf Forest 

Conifer 

2 

Evergreen Broadleaf Forest 

Conifer 

3 

Deciduous Needleleaf Forest 

Deciduous 

4 

Deciduous Broadleaf Forest 

Deciduous 

5 

Mixed Forest 

1/2 Conifer + 1/2 Deciduous 

6 

Closed Shrublands 

1/4 Quartz sand + 3/8 Conifer + 
3/8 Deciduous 

7 

Open Shrubland 

3/4 Quartz sand + 1/8 Conifer + 
1/8 Deciduous 

8 

Woody Savannas 

Grass 

9 

Savannas 

Grass 

10 

Grasslands 

Grass 

11 

Permanent Wetlands 

1/2 Grass + 1/2 Seawater 

12 

Croplands 

Grass 

13 

Urban 

Black Body 

14 

Cropland/Mosaic 

1/2 Grass + 1/4 Conifer + 
1/4 Deciduous 

15 

Snow and Ice 

Mean Of Fine, Medium, and 
Coarse snow and Ice 

16 

Barren 

Quartz sand 

17 

Water Bodies 

Seawater 

18 

Tundra 

Frost 


The decisions on how to associate the surface types were based on the information available, 
and on the authors’ best judgment of how to best characterize the surface type. The 
evergreen needleleaf and the evergreen broadleaf were both assigned the emissivities of the 
conifer sample because no other evergreens were in the archive. The emissivity of the 
deciduous leaf was used for both the deciduous needleleaf and deciduous broadleaf. When 
measurements of other types of trees are available, these emissivity assignments may change, 
although the measurements currently available show that there is little difference between the 
emissivities of evergreen and deciduous forests. Because there is a lack of information on the 
urban surface, a blackbody emissivity of unity was assumed for all spectral regions for the 
urban surface type. The variations in emissivities of dry bare soil is greatest of all surface 
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Figure 3. Surface Types on 10' grid. 








materials and therefore the most difficult to characterize. We were constrained in this version 
of the emissivity maps to use one emissivity classification for all the surface classified as 
barren. Therefore, the barren land was assigned the emissivities of a desert sample composed 
of mostly quartz sand. There have been few measurements made of emissivity of tundra. 
Though Rees (1993) has measured the thermal infrared (8-14 pm) emissivity of a number of 
land cover types in the Svalbard archipelago north of Norway. The observed emissivity 
values fell between 0.941 for sandstone and 0.995 for snow. The window emissivity of frost 
at 0.9806 lies within this range. Thus, the tundra surface type was assigned the emissivity of 
frost. When more is known about the composition and condition of tundra, the emissivity can 
be refined. Surface type 15 is snow and/or ice. To assign emissivity to this surface type the 
emissivities of the 3 snow types were averaged to create an '‘average snow” emissivity and 
then that emissivity was averaged with the emissivity of ice to create the ice/snow emissivity. 
Because the emissivity of ice is very close to 1, the resulting ice/snow emissivity is also close 
to 1. 

The values of emissivity for the 12 Fu-Liou bands, the CERES window and the broadband are 
shown in the Table in Appendix B. Figure 4 shows the band-average emissivities of the 18 
surface types. Figure 5 is the global map of broadband emissivity on the lO' grid. Figure 6 
is the global map of the CERES window channel emissivity on the 10' grid. As additional 
measurements of reflectance and emittance become available, the surface emissivity maps will 
be updated. 

Figure 4. Band-averaged Emissivities for the 18 Surface Types. 
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Figure 4 concluded 
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Figure 5. Broadband Surface Emissivity on 10' grid 
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Figure 6 . CERES Window (8-12 |im) Surface Emissivity on 10' grid 





Computations using emissivity map 


The effect of using the current emissivity maps on surface net longwave fluxes was examined 
using the Gupta longwave model (Gupta et al. 1992). This model is based on a parameterized 
radiative transfer algorithm which computes downward and net LW fluxes at the surface and 
has been extensively validated. This model is being used by CERES for computing surface 
LW fluxes. The model uses meteorology to calculate downward flux. Upward flux from the 
surface is calculated as eoT . The meteorological inputs for the present computation were 
taken from ISCCP-D1 data available on a 2.5° equal-area grid. To run the model it was first 
necessary to change the resolution of the broadband emissivity map, shown in Figure 5, from 
the 10' to a 2.5° equal-area grid. The emissivities from the 1C Y grid were averaged into a 2.5° 
equal-area grid. The resulting broadband emissivity map is shown in Figure 7. The model 
was then run assuming a constant surface emissivity of unity and again with the surface 
emissivity map derived in this work. The resulting net longwave fluxes are shown in Figure 8, 
and the difference in net longwave flux between the two different runs is shown in Figure 9. 
The largest differences (up to 6 Wm ) occur over areas of the Sahara Desert and the Arabian 
Peninsula classified as barren, and open shrubland in Australia. Differences of greater than 
3 Wm" are found over the open shrubland areas of the Western US and Eurasia. There are 
differences of more than 1.5 Wm over the barren areas of Siberia. Because this version of 
the emissivity maps was constrained to have emissivities based only on surface type, the 
emissivity assigned to barren ground was that of quartz sand which is appropriate for the 
equatorial and mid-latitude deserts but may not be applicable to the barren ground in Siberia. 
This treatment of barren ground will be modified as more information becomes available. 

This study has shown that changing the broadband surface emissivity from a constant of 
unity to a variable dependent on surface type can result in differences up to 6 Wm’ in the 
surface net longwave flux. Further refinements and improvements to the emissivity maps are 
necessary to account for the difference in emissivities of barren ground of the quartz deserts 
and the barren ground in high latitudes. 

Future work 

Sand and soil are the surfaces for which it is most difficult to estimate surface emissivities. 
They are also the surfaces for which the emissivities differ most from unity. There is large 
variability in emissivity dependent on composition and surface properties. Emissivity is 
dependent on particle size and soil moisture (Salisbury and D’Aria 1992b). In the case of 
sand there is also a change in emissivity with viewing angle (Snyder et al. 1997). All of the 
aforementioned influences on soil emissivity should be taken into account. The first step in 
improving soil emissivity estimation is to allow for more types of bare soil. The current 
emissivity map assumes quartz sand for the barren surface type because the largest areas of 
barren ground are quartz sand. Plans for modification of the emissivity map are to use the 
Zobler World Soil map (Zobler 1986) in conjunction with the IGBP surface type map. 
Currently there are no laboratory measurments made at wavelengths greater than 16 pm. It is 
desirable to have such measurements in the 18.5-25.0 pm for use with the “polar window” 
band. 

The current emissivity maps are considered to be the basis upon which more refined and 
improved maps will be constructed as information becomes available. When emissivity 
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measurements from MODIS and ASTER are available, they will be incorporated into the 
emissivity maps. Other refinements such as incorporating seasonal surface type and 
vegetation variability, can be made to the emissivity maps before these satellite data are 
available. A surface type can change from green vegetation to brown vegetation and then to 
bare soil over the course of a year. The current map will be modified to take seasonal 
variations into account. There are additional effects on surface emissivity that have yet to be 
considered; however, the current map is a first step in defining more accurate surface 
emissivities. 

Accessing the emissivity maps 


The broadband, CERES window and Fu-Liou band surface emissivity maps may be viewed 
and the data downloaded from the web site: 

http://tanalo.larc.nasa. gov: 8080/surf_htmls/SARB_surf.html 

Sample pages from this web site are given in Appendix C.. 
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Figure 7. Broadband Longwave Surface Emissivity, 2.5 equal-area grid. 
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Figure 9: Difference in NLW Flux (Win' ) : Flux with variable variable e - flux with constant e 





Appendix A. 


The following description of IGBP surface types are adopted from Belward and Loveland 
(1996). 

1. Evergreen Needleleaf Forests: Surface is dominated by trees with a canopy cover of over 
60% and height exceeding 2 meters. Almost all trees remain green all year. Canopy is 
never without green foliage. 

2. Evergreen Broadleaf Forests: Surface is dominated by trees with a canopy cover of over 
60% and height exceeding 2 meters. Almost all trees remain green all year. Canopy is 
never without green foliage. 

3. Deciduous Needleleaf Forests: Surface is dominated by trees with a canopy cover of over 
60% and height exceeding 2 meters. Consists of seasonal needleleaf trees with an annual 
cycle of leaf-on and leaf-off periods. 

4. Deciduous Broadleaf Forests: Surface is dominated by trees with a canopy cover of over 
60% and height exceeding 2 meters. Consists of seasonal broadleaf trees with an annual 
cycle of leaf-on and leaf-off periods. 

5. Mixed Forests: Surface is dominated by trees with a canopy cover of over 60% and 
height exceeding 2 meters. Consists of tree communities with interspersed mixtures or 
mosaics of the other four forest cover types. None of the forest types exceeds 60% of 
the landscape. 

6. Closed Shrublands: Surface consists of woody vegetation less than 2 meters tall and with 
shrub canopy cover of over 60%. The shrub foliage can be either evergreen or 
deciduous. 

7. Open Shrublands: Surface consists of woody vegetation less than 2 meters tall and with 
shrub canopy cover between 10-60%. The shrub foliage can be either evergreen or 
deciduous. 

8. Woody Savannahs: Surface consists of herbaceous and other understory systems, and 
with forest canopy cover between 30-60%. The forest cover height exceeds 2 meters. 

9. Savannahs: Surface consists of herbaceous and other understory systems, and with forest 
canopy cover between 10-30%. The forest cover height exceeds 2 meters. 

10. Grasslands: Surface consists of herbaceous types of cover. Tree and shrub cover is less 
than 10%. 

11. Permanent Wetlands: Surface consists of a permanent mixture of water and herbaceous 
or woody vegetation that cover extensive areas. The vegetation can be present in either 
salt, brackish, or fresh water. 
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12. Croplands: Surface is covered with temporary crops followed by harvest and a bare soil 
period (e.g., single and multiple cropping systems.) Note that perennial woody crops will 
be classified as the appropriate forest or shrub land cover type. 

13. Urban and Built-Up: Surface is covered by buildings and other man-made structures. 
Note that this class will not be mapped from the AVHRR imagery but will be developed 
from the populated places layer that is part of the Digital Chart of the World (Danko, 
1992). 

14. Cropland/Natural Vegetation Mosaics: Surface consists of a mosaic of croplands, forest, 
shrublands, and grasslands in which no one component comprises more than 60% of the 
landscape. 

15. Snow and Ice: Surface is under snow and/or ice cover throughout the year. 

16. Barren: Surface is made up of exposed soil, sand, rocks, or snow which never have more 
than 10% vegetated cover during any time of the year. 

17. Water Bodies: Oceans, seas, lakes, reservoirs, and rivers. The water bodies can be 
composed of either fresh or salt water. 

18. Tundra: Surface is defined by 1GBP to be Barren but is also identified by the Olson 
vegetation map (Olson et al. 1985), as tundra (Arctic wetlands). 


More information on 1GBP is available from the web sites: http://www.igbp.kva.se 
and: http://www.ngdc.noaa.gov: 80/paleo/igbp-dis/index.html 
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Appendix C. 

Sample pages from the web site: 

http://tanalo.larc.nasa.gov:8080/surfJitmls/SARB_surf.html 
For questions or problems involving the site, please contact : 
d.a.rutan@larc.nasa.gov 


26 






Surface & 
Atmospheric 
Radiation 
Budget 



CERES Surface Properties Home Page 


The SARB working group, part of the Clouds and the Earth’s Radiant Energy System CERES 
mission, will calculate profiles of shortwave and longwave fluxes from the surface to the top 
of the atmosphere. The radiation transfer code which will be used was developed by Qiang Fu 
and Kuo Nan Liou, ( Fu & Liou Model ) For proper results the surface boundary condition 
must be specified as a function of the spectral bands of the model based upon the varying 
scenes that the instrument will be observing. For your perusal we’ve placed a set of images 
( Access Here ) that contain the data which will be the starting points for these lower boundary 
conditions. For a more detailed description of how these surface maps are applied in the 
SARB processing consider reading the Surface Properties Description Page . 


Click here for NEWS and updates to these pages. (Latest update: 05/05/99) 


Current Data 

The various data available are listed in the button box below. The images of all but the snow 
and ice data are interactive so if you desire to see an area more closely, aim your pointer in 
the general area you would like to see in detail. If you are interested in a specific latitiude 
and longitude use the Data by Lat/long button to find out all of aour current surface 
information for that location. 
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CERES/SARB Surface Maps for Download 

This page lists all the available information for easy downloading of the maps. All maps, except the digital 
elevation map are 8-bit binary data made on a Sun SPARC Workstation. Their size is 2160 points in 
longitude, 1080 points in latitude or 1/6 degree equal angle. All maps begin at the North Pole, Greenwich 
Meridian. To download a map click on the "MAP" icon or the word "Download", using the right mouse 
button if you're using NETSCAPE. 

* Netscape users, use "save link as" under the right mouse button. ® 


Available CERES Surface Data 



Data Description 

Range 

of 

Values 

Related 

Link 

Last 

Update 

^3 

IGBP+1 CERES scene type map (~2.3Mb) 

1 to 18 


Dec. 03, 
1998 

Download 

Map of Surface Albedo(*100) @ 60Deg Solar 

0 to 100 

Discussion 

Dec. 03, 


Zenith Angle. (~2.3Mb) 



1998 

^3 

Map of Broadband Surface 
Emissivity(* 100).(~2.3Mb) 

0 to 100 


Dec. 03, 
1998 

Download 

Map of Window (8-12Micron) Surface 

0 to 100 

Discussion 

Dec. 03, 


Emissivity(* 100).(~2.3Mb) 



1998 

U3I 

Percentage of water in each 10' grid box.(~2.3Mb) 

0 to 100 


Aug. 01, 
1997 

Download 

Digital elevation in each 10' grid box. (Water 

-500 to 

Discussion 

Aug. 01, 


bodies equal -9999.) (~4.6Mb) 

7000 

(meters) 


1997 

Download 

Snow Map for October, 1986.(~2.3Mb) 

0 to 150 

Discussion 

Aug. 01, 



(inches) 


1997 

US 

Ice Map for October, 1986.(~2.3Mb) 

0 to 100 (%) 


Aug. 01, 
1997 

Download 

Data tables that create the maps. 

- 

Discussion 

May 15, 




1998 

IS! 

Spectral emissivities in the 12 Longwavebands of 
the Fu & Liou code. 

- 


May 15, 
1998 





Back to Surface Properties Home Page 
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